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Abstract. We consider a quantum wire double junction system with each wire 
segment described by a spinless Luttinger model, and study theoretically shot noise 
in this system in the sequential tunneling regime. We find that the non-equilibrium 
plasmonic excitations in the central wire segment give rise to qualitatively different 
behavior compared to the case with equilibrium plasmons. In particular, shot noise 
is greatly enhanced by them, and exceeds the Poisson limit. We show that the 
enhancement can be explained by the emergence of several current-carrying processes, 
and that the effect disappears if the channels effectively collapse to one due to, e.g., 
fast plasmon relaxation processes. 
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1. Introduction 

Shot noise in electronic devices occurs due to the discrete nature of the transported 
charge, and manifests itself clearly in non-equilibrium situations PQI2]- Unlike thermal 
(equilibrium) noise, shot noise provides additional information of the system, which is 
not accessible in an average current measurement [UEIIS1- For this reason, shot noise 
has been a subject of a great number of theoretical and experimental works in recent 
years jE|. 

One of the fundamental issues of current interest has been the effect of coherence 
on the shot noise [Zj. Somewhat contrary to intuition, it has been shown that in most 
cases where comparable theories are available, the semiclassical and quantum mechanical 
descriptions give the same result; it seems that often the shot noise is not sensitive to the 
quantum coherence jBJ. Another important issue related to the shot noise is electron- 
electron correlation effects, which is the subject of the present work. Popular examples 
are the shot noise in a single-electron transistor (SET) in the Coulomb blockade (CB) 
regime jSl EI and in superconducting tunnel junctions [UJ EH HI]- The system of our 
interest in this work is a one-dimensional (ID) interacting electron gas described by the 
Luttinger liquid (LL) model j!2j . There has been a great deal of experimental effort to 
demonstrate that certain low-dimensional systems (e.g., carbon nanotubes) behave as 
Luttinger liquids (THl EH EHj • Shot noise in Luttinger liquid has also been the subject of 
many recent works. For example, the nonequilibrium shot noise of the edge states was 
used to measure the fractional charge of the quasiparticles of the fractional quantum 
Hall states [El E2 EH EH| • 

Very recently, Braggio et al [20] investigated the shot noise in a double junction 
system embedded in a Luttinger liquid. They found some features due to the interplay 
between the CB, the LL behavior, and size quantization. Apart from the peculiar 
quantitative behavior of the shot noise at higher bias voltages, its qualitative behavior as 
a function of bias voltage was similar to that in conventional SET devices. They assumed 
equilibrium distribution of charge- density excitations (in the following called plasmons) 
in all electrodes. On the other hand, in a recent work by some of the present authors 
[21], plasmons in the short central electrode in the absence of environment coupling 
were found to exhibit a highly non-equilibrium distribution. Since shot noise is sensitive 
to non-equilibrium fluctuations in the system, one can expect that the non-equilibrium 
plasmons may affect the shot noise significantly. In this work, we show that the non- 
equilibrium plasmons lead to qualitatively different behavior of the shot noise even at 
relatively low bias voltages. We also show that in the limit of fast plasmon relaxation, 
when the plasmon distribution reaches equilibrium between successive tunneling events, 
the results by Braggio and co-workers are reproduced. 
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2. Model 



We consider a LL/LL/LL double junction with each LL representing a segment of 
quantum wire. We refer to the three segments the left (L) and right (R) leads and the 
dot (.D), respectively. The two leads are adiabatically connected to reservoirs (contacts) 
and the dot (the central short wire segment) is coupled to the two leads by tunnel 
junctions. The L and R leads are described by semi-infinite (spinless) Luttinger models 
while the dot is described by a finite (spinless) Luttinger model |2()|I21|. The low-energy 
transport properties through the system are described by the master equation [2H 122] 
for the probability P(N, {n} , t) to find N (excess) electrons and n m plasmons at mode 
m on the dot at time t. The eigenstates of an isolated dot are denoted by \N, {n}). The 
finite length L D of the dot leads to discrete energy spectrum of plasmons {n} and the 
zero-mode energy associated with N. The equidistant plasmon spectrum is characterized 
by the energy spacing e p = nhvp/ L D g, where g the Luttinger parameter and vp the 
Fermi velocity. The zero-mode energy, characterized by the energy scale Ec = £ P /g, 
accounts for the charging effects on the dot. Electrons are transferred between different 
wire segments by tunneling, and the transition rates Tl/r(N ± 1, {n'} ; N, {n}) between 
states \N, {n}) due to tunneling across the two junctions are obtained using Fermi golden 
rule [2j . Following the notations of Ref. [23 > we consider that the LL-parameter g is 
the same for each wire segment, the capacitances across L/R tunneling junctions are 
taken to be the same (Cx = Cr), but the tunneling amplitudes may differ and their 
ratio is denoted by R = |£l| 2 /I^r| 2 - The explicit form of the transition rates for a wire 
with a single branch of spinless fermions at applied voltage V — Vl — Vr is given by 
Tl,r{N, {n} - N', {n'» = 2n\t L/R \ 2 /h 

xj{E D {N', {n 1 }) - E D (N, {n}) T (N f - N)eV L/R ) (1) 
jN+ i + iPd8n',n-i\N, {n})\ 2 , 
where the contribution of the leads is 



1 1 



7(e) 



2?rA r 



a+ 1 ./3e, 



-f3e/2 



r(« + i) 



and the contribution of overlap matrix elements in the central segment is 
KKM|{n}>| 2 



(2) 



m=l 



n(- 

- LJ - » mg 



n 



l 

L D 

«) 



7rA 

Ld 



X 



lltn - 



T \ n m~ n m\ I 

L ,«) ' 



1 

mg 



(3) 



where n 



«) 



and n 



(» 



Here L^z) is a Laguerre 
polynomial, A a short- distance cutoff, and a = (g^ 1 — 1) the appropriate exponent 
for tunneling into the end of a Luttinger liquid. Note that while the lead contribution 
depends only on the energies E^ of the involved states of the dot, the contribution arising 
from the dot depends on the detailed nature of the participating states. Temperature 
enters the rates through /3 = l/(fcsT); since we are interested in shot noise rather than 
thermal noise, we from now on set ksT to a value far below any relevant energy scale 
in the dot. 
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The plasmons in the two leads, being in contact with reservoirs with many low- 
energy excitations, are assumed to be in equilibria separately. However, the plasmons in 
the dot are driven far from the equilibrium distribution [21] by tunneling electrons. This 
deviation is already significant at relatively low bias voltages 2e p < eV < 2E C . Contrast 
to the equilibrium plasmons [2H] that produce qualitatively similar shot noise as in usual 
SET [3], the non-equilibrium distribution strongly affects the SN of the system. 

The coupling of the system to the environment causes even the plasmon distribution 
in the dot to relax towards an equilibrium. The precise form of the relaxation rate, T p , 
depends on the specific relaxation mechanism, but the physical properties of our concern 
are not sensitive to such details. We use a phenomenological model, obeying detailed 
balance, with 



where W p = £ P J2 m m ( n 'm ~~ n m) is the energy difference of the two many-body states 
with {n 1 } and {n} = n 2 , • • • , n m , ■ ■ •) plasmon occupations, and j3 = l/k B T is 
the inverse temperature. While the plasmon relaxation rate in nanoscale structures 
is difficult to estimate, recent computer simulations on carbon nanotubes indicate 
that plasmon excitations in them have life times of the order of a picosecond, much 
longer than those in three-dimensional structures [2B1- The total transition rate is 
hence given by a sum of rates associated with tunneling and plasmon relaxation, 



T(N', {n'} ; N, {n}) = J2 £=LjR T e {N', {n'} ; N, {n}) + 5 N ,, N T p ({n'} , {n}) where i labels 



the tunnel junctions. The only phenomenological parameter is the time scale of the 
relaxation processes in the dot, 7 P ; all other quantities are determined by the geometry 
of the system and the properties of the tunnel junctions. 

The master equation approach neglects phase coherence of electrons tunneling 
across the two junctions (sequential tunneling regime). The effects of the coherence on 
the SN have been hotly debated in recent years. However, for many mesoscopic systems 
(including low-dimensional interacting systems), the master equation approach yields 
the same results as the fully quantum mechanical description of the SN Therefore, 
we expect that our master equation approach yields qualitatively correct results for the 
system at hand. 

We introduce a matrix notation for the transition rates Y with matrix elements 
defined by 



T p ({n'} , {n}) = 7 P 



W p /e p 



(4) 



e pw P - i 



rt(A0 



J {n'},{n} 



r i (N±l,{n'};N, {n}) , 



(5) 






and 



r P (A0 



- { n '}>{ n } 




{n"} 



Shot noise enhancement 



5 



Using this notation, the master equation is 

|TO> = -P|P(«)> 



(8) 



with r = T p + J2e=L r — — T £ J , where \P(t)) is the column vector (not to be 
confused with the "ket" in quantum mechanics) with elements given by {N, {n} \P(t)) = 
P(N, {n},t). 



3. Shot Noise 

We investigate the shot noise power defined by 

/oo 
dre + ^[(I(t + r)I(t))-(I(t)) 2 } (9) 
-oo 

in the steady state (t — > oo). In terms of the tunneling current Il/r across the 
junction L/R, the total current I is conveniently written as I = Yle=L r{C /Ci)Ii where 
C- 1 = C L l + C R \ The correlation functions K(t) = (I e (t + r)I e (t)) (t — > oo) can be 
deduced from the master equation (jHJ); in the matrix notation [SI El we have 



K(r) = e 2 (^Wl 
N,{n} 



+ Q{-r)h exp(+fr)/ £ + S(t)S u > fr+ 



r; 



|P(oo)> , (10) 



where Q(x) is the unit step function, |P(oo)) is the steady-state solution to the master 
equation (jHJ), and Il/r are tunneling current matrices Il/r = =Fe (j^t/R ~~ ^l/rj- 

For bias voltages and temperatures when only a few levels are involved in the 
transport, the Fano factor can be evaluated analytically, but the general case requires 
numerical integration of equation (jHJ. Figure Q shows a typical behavior of the Fano 




Figure 1. Fano factor F = S(0)/2e (I) as a function of bias voltage (a) for 
g = 0.3,0.5,0.7,1, with no plasmon relaxation on the dot, and (b) for g = 0.5 with 
different plasmon relaxation rates (j p ). Here R = 100 (highly asymmetric junctions), 
and T = 0. 
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factor as a function of bias voltage for different values of the interaction parameter g. 
The dips at eV = 2E C = 2e p /g for g = 0.7 and g = 0.5 in figured] (a) are ascribed to 
the charging effect as typically seen in SET devices [SI IB] • This dip structure was also 
observed by Braggio et al [2D] . 

A remarkable feature that contrasts our model with a conventional SET [HI Ej, 
and the model in Ref. [20], is the enhancement of the shot noise in the bias range 
2e p < eV < 2E C \ see figure [T] (a). As the external bias voltage increases over the 
threshold value 2e p , it provides an electron with enough energy to tunnel across the 
double junction through the plasmon excitation level of mode m = 1, and hence 
the current sharply increases [2D [22J- Unlike a usual resonance channel, however, 
this additional channel is very noisy in the sense that the Fano factor exceeds the 
conventional Poisson limit (F = 1). Therefore, the statistics of the tunneling events 
across the double junction through the plasmon excitations are highly non-Poissonian. 

We emphasize that the SN enhancement is due to the non-equilibrium distribution 
of the plasmons on the dot. To see this, we turn on the plasmon relaxation, 7 P ^ 0. 
For sufficiently large relaxation rates, any plasmon excitations induced by electron 
tunneling relax to an equilibrium before a subsequent tunneling event. Therefore, as the 
plasmon relaxation rate increases, the Fano factor should recover the characteristic of 
the conventional SET or the model of Braggio et al [20] . This is clearly seen in figure ^ 
(b), where the Fano factor is shown as a function of eV for relative plasmon relaxation 
rates 7 p /r = 0, 0.1, 1, 10 4 , where T = \t\ 2 /h 2 v F L D with \t\- 2 = \t L \~ 2 + \t R \~ 2 . 

A simple explanation of the enhanced noise may be obtained by considering a 
situation in which only three states of the dot are energetically allowed. We denote the 
states |0), \1A) and \1B) where the integer denotes the (excess) number of charges, and 
A and B label different states with the same particle number^. The transition rates 
between these states, caused by tunneling event into and out of the dot, are given by 
where s = ± indicates if the transition is from charge state to charge state 1 or 
vice versa, and the label a = A, B tells which of the states |L4) and \1B) is involved 
in the transition. Hence, there are two current carrying processes involving transitions 
|0) — ► | la) — ► |0), and, in a time sequence, the two processes alternate randomly. This 
extra degree of freedom results in additional noise, and is reflected in the Fano factor 
(which can be calculated as described above) 



that can exceed unity if the two processes correspond to very different current levels; in 
the case studied above this happens when the interaction is sufficiently strong. If there 
is a fast internal relaxation process that allows the state 1 to decay to |L4), the two 

| In the case studied above, this three-state model is reasonably accurate for 2e p < eV <C 2Ec, when 
A corresponds to the state with no plasmon excitations and B to the state with one plasmon at the 
mode m = 1 (ni = 1). 
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Figure 2. Finite frequency shot noise S(w)/2el as a function of the bias eV/Ec and 
frequency eoj/I for R = 100 and g = 0.5, at Nq = 1/2 (T = 0), (a) with no plasmon 
relaxation (j p = 0) and (b) with fast plasmon relaxation (j p = 10 4 ). In (c), the longest 
relaxation time with 7 P = as a function of the applied voltage for strong and weak 
interactions (g = 0.5 and g = 0.7, respectively). The relaxation time is given in units 
of T = \t\ 2 /h 2 v F L D with \t\- 2 = \t L \- 2 + \t R \- 2 . 



channels effectively collapse into one, and the Fano factor is given by 

F - K+n+r+P e| 2' 11 (12) 

which is the result for a single channel system with transition rates T~ = and 
T + = + T^. Hence, the presence of charge-conserving internal relaxation processes 
reduces the effective number of transport processes, and the noise associated with the 
random alternation between them. 

The shot noise enhancement persists to quite high frequencies as shown in figure 121 
where the shot noise power S(u) / (2el) is shown as a function of the bias and frequency 
for a strong interaction (g = 1/2) and asymmetric tunnel barriers (R = 100) without 
plasmon relaxation (a) and with fast plasmon relaxations (b). 

The frequency dependence of noise reveals two additional consequences of non- 
equilibrium plasmons. Firstly, the asymptotic high-frequency value of S(uj) is increased 
by the non-equilibrium distribution, and secondly, the characteristic frequency scale of 
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the noise is seen to depend quite sensitively on the applied voltage. 

In the high frequency limit, uj — > oo, the correlation effects are lost in the 
noise power except the <5(r)-term in (JTHJ) that reflects the Pauli exclusion |HJ, and the 
asymptotic value of the noise spectrum reduces to 

^ — > = 2e ((0k? + <fe) = > + ^ (13) 

where A L / R in the steady-state is given by 

A e = eJ2 (N, {n}\ft + f;|P(oo)> . (14) 

N,{n} 

The quantity Ai counts the total number of tunneling events across the junction £, 
without regard to direction. In the absence of tunneling events against the voltage, 
which is the case e.g. in the limit of fast plasmon relaxation and low temperatures, 
the limiting value of S(u) is simply proportional to the current. In the non-equilibrium 
case, however, some plasmon modes with sufficiently high energy are occupied so that 
tunneling against the voltage is not negligible, which results in the enhanced high- 
frequency noise above the customary limiting value S(u — ► oo) = el. 

The frequency dependence of S(u) comes from terms of the form XjJ {u 2 + A|) (see 
Equation (jlOj) ). where Xj are the non-zero eigenvalues of the transition matrix T, that 
is, the characteristic relaxation rates of the system. The longest relaxation time, shown 
in figureEfc) for g = 0.5 and g = 0.7, exhibits maxima at threshold voltages when more 
states (either charge states or additional plasmon excitations) become involved in the 
transport. The long relaxation time reflects the slow transition rates to, or from, these 
states near the energy threshold. For strong interactions (g = 0.5 in the figure) clear 
maxima are seen both at thresholds for additional plasmon states (eV = (2g)Ec) and for 
new charge states (eV = 2Ec), while for weaker interactions, the structure at plasmon 
excitations is almost washed away. This is due to the fact that for weak interactions, the 
transition rates are roughly given by step functions, in contrast to the strong interaction 
case when they increase like a power law with a larger exponent. 

4. Conclusions 

We have investigated the shot noise of a Luttinger-liquid double-barrier structure. The 
non-equilibrium plasmons on the central segment affect strongly current fluctuations 
even at low bias voltages 2e p < eV < 2Ec- The shot noise is enhanced beyond 
the Poissonian limit due to the emergence of competing transport processes, the 
enhancement is more pronounced for stronger interactions, and persists to unexpectedly 
high frequencies. 
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